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Abstract. We discuss the spectral energy distribution of three very low mass objects in Chamaeleon I for 
which ground-based spectroscopy and photometry as well as ISO measurements in the mid-infrared are avail- 
able (Comeron et al. 2000; Persi et al. 2000). One of these stars (Cha Hal) is a bona-fide brown dwarf, with 
mass 0.04-0.05 Mq. We show that the observed emission is very well described by models of circumstellar disks 
identical to those associated to T Tauri stars, scaled down to keep the ratio of the disk-to-star mass constant and 
to the appropriate stellar parameters. This result provides a first indication that the formation mechanism of T 
Tauri stars (via core contraction and formation of an accretion disk) extends to objects in the brown dwarf mass 
range. 



1. Introduction 

The last few years have seen an enormous progress in our 
understanding of sub-stellar mass objects, as more and 
more objects in the mass range of brown dwarfs (BD; 
M^,<0.75 Mq) and giant planets (< 0.015 Mq) are found, 
both in the field and in regions of recent star formation 
(Basri 2000; Lucas & Roche [2000 Zapatero Osorio et 
al. 200C| ). It is now clear that "free-floating" very low mass 
objects exist. Do they form, as low-mass stars do, from the 
collapse of a molecular core? This, although apparently 
very reasonable, is not the only possibility. Alternative hy- 
potheses have been suggested; for example that BD form 
in gravitationally unstable regions of protostellar disks 
(Pickett et al. 2000) or that they are stellar embryos. 



whose further growth is prevented by dynamic al ejec tions 
from small stellar systems (Reipurth & Clarke 2001 ). 

A full understanding of the formation mechanism(s) of 
sub-stellar objects will take time. At present, it is timely 
to begin to explore the properties of BD systems in re- 
gions of star formation and compare them to the much 
better known pre-main-sequence low-mass stars, the T 
Tauri stars (TTS). The most important clue to a forma- 
tion mechanism involving accretion from a parental core is 
the presence of a circumstellar disk. Claims of evidence of 



excess emission in the near (Oasa et al. 1999; Muencli et 
al. pOOlD and mid-IR (Comeron et al. pOOOj ) in BD or BD- 
candidates in star-forming regions are beginning to appear 
in the literature. Hq, (Comeron et al. 1999; Comeron et al. 



200C) and X-ray emission (Neuhauser & Comeron 1998; 
Comeron et al. 2000] ) are seen in some very low mass 
objects. MuzeroUe et al. (pOOO) detect evidence of mag- 



netosphcric accretion at a very low rate '--^ 5 x 10^^^ Mq 
yr^^ in the spectrum of the M6 object V410 Anon 13, 
whose mass is estimated in the range 0.04-0.06 Mq. 

In very cold BD, the most convincing evidence of 
excess emission, the accepted signature of circumstellar 
disks, can only be obtained in the mid-IR. There are 
at present very few bona-fide BD (i.e., stars with spec- 
troscopical classification) with measured mid-IR excess. 
Comeron et al. (2000) list a small group of BD detected 
in the mid-IR survey of Cham I at 6.7 and 14.3 fim by ISO 
(Persi et al. 200C). Of these, only three have detections in 
both bands, the one with the lowest mass being Cha Hal, 
a M7.5 star with mass 0.04-0.05 Mq. 

In this paper, we show that disk models analogous to 
those developed for TTS can well account for the observed 
spectral energy distribution (SED) of these three stars, 
providing strong support to the idea that BD form like 
the more massive TTS. We will also speculate somewhat 
on the derived disk properties, and on the possibility of 
extending this kind of study to objects of lower mass. 
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2. Disk Models 

We have computed the emission expected from a circum- 
stellar disk heated by the irradiation of the central star fol- 
lowing the method outlined by Chiang & Goldreich ( |1997| 
CG97), which has been successfully apphed to pre-main- 
sequence TTS and Hcrb ig Ac stars (Natta et al. [2000a| 
^OOll ; Chiang et al. |200l|) . CG97 consider flared disks, in 
hydrostatic equilibrium in the vertical direction, and make 
a number of simplifying assumptions which permit the 
computation of the resulting SED in a quick and efficient 
fashion. Although not entirely self-consistent, such models 
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provide a good first approximation to the SED, more than 
sufficient for the purpose of this paper. 

We have taken for the various model parameters val- 
ues typical of pre-main-sequence stars, scaled down where 
necessary. We have assumed that the circumstellar disk 
extends inwards to i?*, has outer radius Ru= 100 AU, 
total mass Md=0.03 M^(Natta et al. ^000h| ), power-law 
surface density S ex i?~^'^; the dust in the disk midplane 
ha s opa city k — 0.01(1.3mm/A) cm^ g~^ (Beckwith et 
al. I990I ). On the disk surface, we assume the mixture of 



carbonaceous and silicate grains that provides a good fit 
to the mid-I R em ission of pre-main-sequence stars (see 
Natta et al. ^001 for details). Most of these parameters 
are either irrelevant for the determination of the mid-IR 
flux, or appear in combinations, and cannot be individu- 
ally constrained by the d ata available (see, for a discus- 
sion, Chiang et al. 2001 ). At this stage, only the most 
"standard" assumptions are justified. 

The stellar properties have been determined by 
Comeron et al. ( |200C ). The most important parameter 
for the disk SED calculation is the stellar luminosity 
and, to a lesser degree, the ratio M^^/Ri,, that controls 
the disk flaring angle (roughly oc(i?*/M*)4/7; CG97). We 
have used in displaying the results of our calculations the 
model stellar atmospheres of AUard et al. ( ^000| ; 2001). 
When comparing them with the broad-band visual and 
near-infrared photometry, we found good agreement only 
for values of the effective temperatures significantly lower 
than the spectroscopically determined values of Comeron 



et al. (200C). We do not know if this effect has any signif- 
icance, given the uncertainties on the temperature scale 
in this mass range. In any case, our effective temperatures 
are well within the range of values expected for field brown 
dwarf s of s imilar spectral type (see, for example, Leggett 
et al. 2001), and their exact value is unimportant for the 
disk SED determination. The extinction and luminosity 
we derive are identical to the Comeron et al. values for 
Cha Hal and Cha IIq;2, while they are definitely higher 
for Cha IIa9, for which we estimate a luminosity of 0.018 
1/0 as compared with Comeron et al. value of 0.0056 Lq. 
The adopted values of the stellar parameters are given in 
the figure captions. We have not re-determined the stellar 
masses, for which we used the Comeron et al. estimates. 

The most interesting of the three stars is Cha Hal 
which has the lowest mass of the three (about 0.04-0.05 
Mq). The results are shown in Fig. H. The disk predic- 
tions fit extremely well the observed points at all available 
wavelengths. 

The results for the other two stars are shown in Fig. 
and ^, respectively. Our standard disk model fits rather 
well the observed photometric points, although we note 
a tendency of the 14.3 /itm point to lie slightly below the 
model predictions. This discrepancy, however, is only 20% 
for Cha Ha2 and 25% (2cr) for Cha Ha9. For this star, 
we show in Fig. 3 the predictions of the same disk model 
seen with an inclination angle of 75 deg (dashed line). 
The flux is only slightly reduced, since it is dominated at 
most wavelengths by emission of optically thin material 
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Fig. 1. SED of Cha Hal. The star has T,=2400 K, 
L*=0.01 Lq, M^,=0.05 Mq, extinction in the I band 
Ai=0.11 mag. Disk parameters as in the text, the viewing 
angle is 0° (face-on). The jagged fine shows the AUard et 
al. (2001) model atmosphere, the solid fine the disk model 
predictions. The short and long dashed lines show sep- 
arately the contribution to the SED of the disk surface 
and midplane. The filled squares with error bars are the 
observed fluxes from Comeron et al. (2000). 



(see Fig. || and g). The dot-dashed line shows the SED 
predicted by a geometrically flat disk seen face-on (dot- 
dashed line). While the difference at long wavelengths is 
very large, the mid-infrared flux would still be consistent 
with the ISO measurements. One point worth to notice is 
that in these low luminosity objects the mir-IR is emit- 
ted by the inner disk, and the ISO points would not be 
consistent with a disk inner hole larger than about 3 -R*. 
The results shown in Fig. 1, 2 and 3 provide good ev- 
idence that the same kind of disks which reproduce the 
properties of pre-main-sequence stars exist around lower- 
mass objects, including a bona-fide BD such as Cha Hal. 
In all three cases, the disks need to be optically thick in 
the mid-infrared. This, however, sets only a weak con- 
straint on the mass of the disk, which remains optically 
thick at 14.3 /im as long as Md>1.5 x 10^®Ki4.3^m Mq, 



where k 



14.3/im 



is in units of cm^ g ^. For K14 



3fim 



3-10 



(Henning & Stognicnko 2000), this limit translates into 
Md>10~^ Mq, or >2 x 10""* A/*, about ten times lower 
than the smallest measur ed ratios in pre-main-sequence 
stars (Natta et al. 2000b). The disk mass, and, to some 
extent, its size, can be better determined from millime- 
ter data. For a typical value A/d~0.03 Mi,, we predict for 
Cha Hal a 1.3mm flux of 3 mJy, well within the range of 
existing millimeter telescopes. 
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Fig. 2. Same as Fig. 1 for Cha Ha2. The star has ^,=2550 
K, L*=0.035 Lq, M^=0M Mq, Ai^OA mag. The disk 
viewing angle is 0°. Note that the formal errors on the mid- 
infrared fluxes quoted by Comeron et al. are 6% and 5% at 
6.7 and 14.3 /zm, respectively, while we have plotted error 
bars corresponding to a more realistic 10% uncertainties. 
The point at 3.6 /im is from Kenyon & Gomez 2000. 
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Fig. 3. Same as Fig. 1 for Cha Ha9. The star has r^=2500 
K, L^=0.018 Lq, M*=0.09 Mq, Ai=1.5 mag. The disk 
viewing angle (solid line) is 0°. Also shown are the pre- 
dictions of the same disk seen at an inclination of 75° 
(dashed line) and of a flat disk seen at 0° (dot-dashed 
line). To avoid confusion, we do not show the disk surface 
and midplane constributions separately. 



Cha Hal has spectral type M7.5 and mass of 0.04-0.05 
Mq. Of course, one would hke to extend our knowledge of 
circumstellar disks to even lower mass objects, such as are 
currently being discovered in the Orion Nebula Cluster 
(ONC; Lucas fc R oche 2000|) and in a Orionis (Zapatero 
Osorio et al. |2000| ). Fig. 4 shows the predicted SED for 
objects of lower and lower mass, as labelled, assuming an 
age of 2 Myr. Photosphcric effective temperatures and lu- 
minosities are from Baraffe et al. (1998) and Burrows et 
al. ( 2001 ). Given the heuristic purpose of this figure, we 
have assumed that the stellar emission can be described by 
a black body at T*. This, as already noted, does not affect 
the calculations of the disk emission, but is a poor de- 
scription of the photospheric spectrum. As a consequence, 
the SEDs in Fig. 4 are not realistic at wavelengths shorter 
than ^4 /im, where the photospheric emission dominates. 

Our calculations show that in nearby star forming re- 
gions it will be possible to detect disks around young sub- 
stellar objects with current instrumentation. Modern mid- 
infrared cameras at large telescopes have a 10 /xm sen- 
sitivity allowing detection of disk emission from systems 
more massive than '-~^10 Mj. The upcoming space missions 
(SIRTF, HERSCHEL and especially NGST) will allow to 
detect disk emission around planetary-mass objects in the 
ONC and the a Orionis clusters. Dust spectroscopy in the 
mid-infrared, as obtained by ISO for luminous pre-main- 
sequence A stars (see Waelkens et al. 199(]| ) , is beyond the 
capability of even the largest ground-based telescopes, but 



for the most massive and nearby BD; in principle SIRTF 
and NGST will allow the observation of the spectra of 
lower mass or more distant objects, although the former 
satellite may be seriously affected by confusion problems. 
Perhaps the most interesting observations, determining 
the amount of circumstellar material around these objects 
and possibly probing the disk kinematics, will be those in 
the millimeter-wave range. As noted above, our disk model 
for Cha Hal predicts a continuum flux (about 3 mJy at 
1.3 mm) well within the range of detectability with cur- 
rent instrumentation. To detect disks around less massive 
and more distant systems and to attempt the detection of 
molecular line emission in order to study the disk kine- 
matics, we will need to wait for the ALMA array to be 
operational. 

3. Conclusions 

The presence and properties of disks around very low mass 
objects, which are being discovered in regions of star for- 
mation, are of crucial importance for understanding their 
formation. If disks exist, and if their properties are analo- 
gous to those of disks around TTS, we can conclude that 
the formation mechanism of TTS, by collapse and accre- 
tion of a molecular core, is very likely to extend to the 
lowest mass objects. 

This paper presents the first attempt to test quanti- 
tatively this hypothesis. We have modeled the emission 
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Fig. 4. Plot of the expected flux for BD of decreasing 
mass, from 40 Jupiter mass (0.04 Mq) to 5 Mj, as labelled. 
Below ~4 /xm, the emission is dominated by the stellar 
photosphere, assumed here as a black body at T^. The 
expected detection limit (lOcr in 1 hr) are shown for the 
VLT/VISIR mid-infrared camera and for the far-infrared 
bolometers on HERSCHEL (horizontal bars). 

expected from disk models of properties identical to those 
of TTS for three objects in Chamaeleon I. These objects 
have been studied spectroscopically by Comeron et al. 
(2000), who attributed to them spectral types M7.5-M6 
and masses in the range 0.04-0.09 Mq. The lowest mass 
object (Cha Hal) is certainly a bona-fide BD. These ob- 
jects are of particular importance because they have been 
detected at 6.7 and 14.3 /xm by ISO (Persi et al. 2000), 
providing a unique (until now) possibility to constrain the 
disk properties of BD. Inspection of Fig.l to 3 shows that 
it is extremely difhcult to infer the existence of disks from 
near-infrared photometry, where the emission is largely 
dominated by the stellar photosphere. 

In all three cases, and in particular in the most inter- 
esting object Cha Hal, optically thick, flared disks are 
required to account for the mid-infrared fluxes. The disks 
are heated by the central star, and we expect that rela- 
tively small grains on the disk surface contribute signif- 
icantly to the observed mid-infrared emission. If higher 
spectral resolution data in this range were available, we 
predict that one would observe the 10 /zm silicate feature 
in emission. Such disks are identical to those around TTS, 
just scaled to the appropriate stellar parameters. 

This result seems to us very interesting. It provides a 
first indication that the TTS formation mechanism (from 
a collapsing core via an accretion disk) extends to objects 
of 0.04-0.05 Mq. It does not, however, rule out the pos- 
sibility that BD are ejected stellar embryos, as suggested 



by Reipurth & Clarke (2001). In their model, the embryos 
may keep a small circumstellar disk, of few AU size. The 
existing data, limited to wavelengths shorter than '^15 
yum, can only set a limit to the disk radius of i?D^l AU. 
There are two ways to proceed, the first and more obvious 
is to detect the millimeter emission of at least some BD, 
which will be well below detection in the embryo hypoth- 
esis. The second is to search for other BD and planetary 
mass objects in the mid-infrared. A high frequency of disks 
is not predicted by the Reipurth & Clarke model, since the 
small, truncated disks, associated with the embryos, not 
fed by any surrounding core, will rapidly disappear. 
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